Vortex-vacancy interactions in two-dimensional easy-plane magnets 
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For a model with isotropic nearest neighbor exchange combined with easy-plane exchange or 
single-ion anisotropics, the static effects of a magnetic vacancy site on a nearby magnetic vortex are 
analyzed on square, hexagonal and triangular lattices. Numerical energy minimization and linear 
stability analysis using the vortex instability mode are employed. When the vortex is centered on 
a vacancy, the critical anisotropies where the stable vortex structure switches from out-of-plane to 
planar form are determined, and the vortex energies and magnetizations are found as functions of 
anisotropy. Consistent with square lattice calculations by Zaspel et al, the strength of anisotropy 
needed to stabilize a vortex in the planar form is reduced when the vortex is centered on a vacancy, 
for all three lattices studied. The vortex-on-vacancy energy is found to be smaller than the typical 
energy of a vortex centered between lattice sites in a system without vacancies. For a vortex 
separated from a vacancy, the energy found as a function of separation demonstrates an attractive 
potential between the two. 

PACS numbers: 75.10.Hk, 75.30.Ds, 75.40.Gb, 75.40.Mg 
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I. INTRODUCTION 

Magnetic vortices in two-dimensional (2D) systems 
with easy-plane coupling have attracted a number of 
investigations related to the structure and stability of 
the in-plane and out-of-plane vortex types liiSi^iiiSiSii The 
presence of these two vortex types- was found in nu- 
merical simulations to exhibit an instability at a criti- 
cal anisotropyfSi^ at which one type can transform to the 
other due to favorable energetics. The instability has 
been explained in terms of the energy balance and soft- 
ening of a dynamic mode on a vortes^ using discrete lat- 
tice models, with precisely located critical points found 
by analyzing progressively larger and larger vortex core 
regions&i 

A magnetic material is likely to have some distribu- 
tion of magnetically inert or nonmagnetic impurities at 
some of the lattice sites. Effectively these are magnetic 
vacancies with respect to the magnetic degrees of free- 
dom. Each vacancy site is expected to play a role in 
the structure and dynamics of all the magnetic excita- 
tions present. Zaspel et al^ considered a model for such 
a nonmagnetic impurity in a square lattice 2D ferromag- 
net (FM). Applying the same discrete core methods as 
was done for uniform systems, they studied the structure 
and stability of a vortex centered on a vacancy site on a 
square lattice. They found that the easy-plane anisotropy 
strength necessary to keep the vortex in the in-plane form 
was drastically reduced by the vacancy. The effect is sig- 
nificant: the critical anisotropy strength 5c (as a fraction 
of the exchange energy, see definition and Hamiltonian 
below) was found to change from about 0.2966 for a free 
vortex to a new value 0.0429 when pinned on a vacancy. 
It means that in real materials even the presence of a 
low density of vacancies could dramatically influence the 
static and dynamic magnetic properties. Further studies 
have indicated significant dynamic features attributed to 
vortices pinned on impurities in isotropic antiferromag- 



nets (AFM),2iiaiii 

On the other hand, M61 et al"'^^ made a continuum cal- 
culation for a vortex near a nonmagnetic impurity un- 
der the assumption that the impurity produces a cer- 
tain global deformation of the vortex spin structure. The 
calculation indicated an effective repulsive potential be- 
tween a planar vortex and the nonmagnetic impurity. 
This would suggest that a spin vacancy could not at- 
tract or pin a nearby vortex. Then the discrete core 
model for the calculation of critical anisotropy would 
be of questionable value. Monte Carlo simulations of 
a planar rotator model-'^^ indicated a lowering of the 
Berezinskii-Kosterlitz-Thouless transition temperature- 
due to the presence of magnetic vacancies. Using a con- 
tinuum model, Leonel et aU'^ concluded that a nonmag- 
netic impurity will repel either an individual vortex or 
even a vortex-antivortex pair, and that this effect can 
even force the BKT transition temperature to zero at a 
critical density of vacancies. On the other hand, study of 
a 2D isotropic Heisenberg antiferromagnet by M61 et ali^ 
determined oscillatory dynamic modes of solitons pinned 
to vacancies, which requires an attractive restoring po- 
tential. These are contradictory results, therefore, fur- 
ther investigation of the vacancy-vortex interaction po- 
tential is needed. 

Although these are related but distinct models (three- 
component easy-plane, one-component planar, and AFM 
vs. FM), it would be satisfying to understand better the 
vortex-vacancy interaction problem in general, partly to 
resolve this discrepancy, and to gain a better understand- 
ing of the overall energetics and the influence of different 
lattices. The studies here will involve only lattice mod- 
els; they may give some insight when compared to the 
above mentioned continuum predictions. Ultimately a 
lattice model for a vacancy is the original source of any 
continuum theoretical description; the continuum theory 
cannot completely describe subtle differences ocurring on 
the different lattices at short length scales. In this re- 
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gard, the vortex-vacancy interaction is studied here on 
2D square, hexagonal and triangular lattices. 

For simplicity, we consider the effect of a single va- 
cancy on a single magnetic vortex. The static struc- 
ture and properties of the vortex are found when the 
vortex is centered on the vacancy. Additionally, an ef- 
fective vortex-vacancy interaction potential is estimated 
as a function of their separation. It is found that the 
vortex-vacancy interaction potential is always attractive, 
and that at well-chosen anisotropy, the dynamic process 
of pinning of a vortex onto a vacancy could lead to strong 
out-of-plane spin fluctuations at very low energy cost. 



II. EASY-PLANE MODEL 

The model to be investigated has classical spins defined 
at lattice sites n interacting with nearest neighbors at 
displacements a according to the Hamiltonian, 

H = J 2 ^ ['^n'^n+a+'^'n'S'n+a+'^'^'n'S'n+a +'^('5'n)^} 

n a 

(1) 

Parameter A determines the exchange anisotropy 
strength 5 = 1 — \ and d is the single ion anisotropy 
strength. Positive values of 5 and d correspond to the 
easy-plane anisotropy that is of interest here. The z- 
axis is the hard axis; xy is the easy plane. A ferro- 
magnetic interaction is assumed. For static vortex prop- 
erties, there are essentially no differences in energy or 
critical anisotropics for the antiferromagnetic model on 
square and hexagonal (bipartite) lattices. The vortex 
structure in terms of staggered magnetization for AFM 
vortices is the same as that for FM vortices with the 
same anisotropy parameters. Obviously these comments 
do not apply to the AFM model on the triangular lattice, 
due to each site having six nearest neighbors, requiring 
three sublattices, which results in frustration. With this 
exception, study of the FM model is sufficient for deter- 
mining the static vortex properties. Dynamic properties, 
however, will be different for AFM vs. FM models. 



angles lose the perfect circular symmetry of this for- 
mula, and obtain modifications largest near the vortex 
core. They satisfy a discrete lattice nonlinear Laplace 
equation. 
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(The same equation will give the vortex structure also 
when a vacancy exists.) The out-of-plane component m 
satisfies an appropriate energy minimization equation re- 
sulting from H . 

It is well known that two vortex types (out-of-plane, 
at weak anisotropy, and in-plane, at stronger anisotropy) 
can be present in this model. The out-of-plane vortex 
has a nonzero profile for m„ whose magnitude typically 
peaks somewhere near the vortex core, and diminishes 
over a length scale determined by the inverse anisotropy 
strength. The in-plane or planar vortex has m„ = 
everywhere. The in-plane vortices require a minimum 
anisotropy strength (denoted by critical parameters Ac 
or dc for the different forms of anisotropy) in order to 
be stabilized. Conversely, the out-of-plane vortices are 
stable only if the anisotropy is weaker than the critical 
value. For a vortex centered in a plaquette (between 
lattice sites), the critical anisotropy strength for in-plane 
vortex stabilization increases with the coordination num- 
ber (z) of the lattice.® For hexagonal {z = 3), square 
(z = 4) and triangular (z — 6) lattices, the respec- 
tive critical values are 5c ~ 0.1670,0.2966,0.3871 and 
dc w 0.2212, 0.4690, 0.8391. Additionally, there is a weak 
attractive potential pulling the vortex center to a posi- 
tion of high symmetry at the center of a plaquette. In 
the present study, we investigate how the presence of a 
vacancy could modify the critical anisotropics for the dif- 
ferent lattices, and how a vacancy otherwise affects the 
vortex structure and energetics. 



III. VORTEX ON A VACANCY 



A. Vortices in a Uniform System 

It is convenient to write the spins of length S in terms 
of in-plane angles 0n = tan~^(5-(/S'^^) and scaled out-of- 
plane component to„ = S^/S, 

Sn = S{y/l - TO^ COS0n, y/ 1 - siu (/)„ , TOn ) . (2) 

In this notation a static vortex at position fy = (xy, yy) in 
continuum theory has in-plane angle in the configuration 

0(r) = ^^{r) = tan-i f ^3^) + ^0, (3) 
y X Xy J 

where (/)o is an arbitrary constant, set to zero here with- 
out loss of generality. On a lattice, the in-plane vortex 



A system is considered with a single magnetic vacancy 
at its center, which is taken as the origin of coordinates, 
(0,0). The center of the vortex is coincident with the 
vacant site. For the purpose of analysis and numerical 
calculations, a finite circular system of radius i? is con- 
sidered, on either hexagonal (i.e., honeycomb), square or 
triangular lattices. This means lattice sites are set up 
surrounding the origin, according to the chosen lattice, 
and only those within radius R are kept. The system is 
abruptly cut off at radius r — R] effectively it is a sys- 
tem with a free boundary condition at this radius. The 
missing magnetic site corresponds to a number of missing 
exchange bonds equal to the coordination number of the 
lattice, z = 3, 4, 6 for hexagonal, square or triangular lat- 
tices. Thus, the order of magnitude of the energy change 
due to a vacancy should be proportional to zJS"^ . 
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A. Critical Anisotropy Calculations 

The critical anisotropics are determined by locating 
the values where it becomes energetically favorable for 
the in-plane vortex to develop nonzero out-of-plane com- 
ponents. The method used in Refs. 0, considering a set of 
core spins that is most important in the energetics, is ap- 
plied here. Each spin is assumed to have small amplitude 
dynamic deviations from the in-plane vortex structure: 

(j>n = <i>l + ifin, = Srrin. (5) 

The linearized equations of motion for the fluctuations 
if and m produce a zero-frequency normal mode when 
the anisotropy strength is at a critical value. The mode 
obtains an imaginary frequency when the anisotropy 
strength is weaker than the critical value, signifying the 
growth of nonzero out-of-plane spin components, and the 
transition to an out-of-plane vortex structure. 



1. Ring Coordinates on Different Lattices 

This special vortex instability mode of oscillation has a 
nearly circular symmetry around the vortex core, which 
allows for great simplification of its analysis. In fact, cir- 
cular symmetry is not exactly valid on a lattice. Instead, 
one supposes that a given site with coordinates (a;n,?/n), 
and other sites related by symmetry operations appropri- 
ate to that lattice, all have the same in-plane and out-of- 
planc spin fluctuations. These are refered to as "rings" in 
Ref. jji. All the sites of a given ring have the same radius, 
however, all sites at the same radius do not necessar- 
ily belong to the same ring. Also, because the vortex is 
assumed to be centered on the vacancy, rather than cen- 
tered in a plaquette, the ring definitions are different for 
the vacancy problem than for the uniform system. 

On a square lattice, the rings can be defined using 
principal sites in the first octant of the xy plane (Fig. 
0. Around the vacant origin, using lattice constant a, 
a site at {x,y) = a(i,j) and all its symmetry related 
sites at o(±z,±j) and a(±j,±i) compose a ring. A ring 
a is defined by a pair of positive integers, a = 
The allowed rings have i = 1,2,3... and j = 0,1,2... 
with j < i; the restriction i ^ produces the vacancy. 
The number of members in a ring will be denoted as /icK ■ 
The rings with j — or j — i lie along lines of high 
symmetry and have only fj.a = 4. All other rings have 
/Xq = 8. Unlike the uniform system problem, there are no 
exchange interactions within an individual ring (no intra- 
ring exchange energy, E^ca = 0). For neighboring rings 
a — and a' = the number of exchange bonds 

between them is denoted by Ca,a' ■ One has Caa' — 4: — z 
when j = j' — (both rings' principal sites on the x-axis, 
a line of high symmetry), and Caa' = 8 = 2z otherwise. 
Before relaxation to a minimum energy in-plane vortex 
state, the in-plane vortex angle for principal ring site with 
is simply = tan'^ [j/i). 



On a hexagonal lattice, a set of principal sites within 
the first sextant of the xy plane define the rings (Fig. 
|2J|. With the vacancy at the origin, and using integers 

1 — 1,2,3... and j — 0,1,2... with j < i, the sites on a 
ring lie at principal sites {x,y) — a{^i, ^j), and other 
sites related by 120° rotations and reflections across the 
line i = j. The restriction that {i + j) is even produces a 
triangular lattice and the additional restriction i — 3j — 

2 7^ 6n, n ^integer, removes the hexagon center sites. 
One can note that the principal ring sites might also be 
constructed as the positions r = la + jb, where I = ^{i — 
j), using basis vectors for hex/tri lattices, a = a(l, 0) and 

b = a(i, -^). The rings along major lines of symmetry 
where j — or j = i have number of members /Xq. = 3; all 
other rings have ii^ = G- There is no intra-ring exchange 
energy, Ecc,a = 0. Neighboring rings a = and a' = 
along major lines of symmetry j ~ j' = (both 
along a axis) or with (i — j) = {i' — f ) (both along b axis) 
have number of exchange bonds Caa' = 3 = z. Other 
neighboring rings have Caa' = 6 = 2z. Before relaxation, 
the in-plane vortex angle for principal ring site is 
$^ = tan"i(\/3jA). 

For the triangular lattice system, a set of principal 
sites in the first 30° above the x-axis is sufficient to define 
the rings (Fig. One needs ring sites as on the hexag- 
onal lattice, with {x,y) = a(^i, -^j), where i — 1,2,3... 
and i = 0, 1, 2... with the restriction that {i + j) is even. 
To constrain to the principal sites in the first 30° slice 
requires j < The rings with principal site along high 
symmetry lines at 0° {j — 0) and 30° {j = ^i) to the 
-|- X-axis have number of members /Iq = 6; all others have 
/iQ = 12. The rings such as (i,j) = (2, 0), (5, 1), (8, 2), 
etc., whose principal sites lie just below the 30° line, sat- 
isfying j = i(z — 2), have six intra-ring exchange bonds 
across the 30° line. The intra-ring exchange energy in one 
of these rings with principal site in-plane vortex angle $^ 
and out-of-plane fluctuation rua is 

£;ee,a = {(1 - ml) cos[2(30° - $„)] + \ml]. 

(6) 

Just as for the hexagonal lattice, neighboring rings a — 
(«,j) and a' = (*'ij') along major lines of symmetry 
i = j' = or with (i—j) = {i'—j') have Caa' = 6 = z, and 
other neighboring rings have Caa' — 12 — 2z. The unre- 
laxed in-plane vortex angle for principal ring site is 
=tan-i(y3j7z). 



2. Relaxed Static In-Plane Vortices 

It was mentioned in the Introduction that the static in- 
plane angles satisfy the discrete nonlinear Laplace equa- 
tion , even with a vacancy in the system. This is seen 
even from minimizing the total energy in the ring coordi- 
nates, under the stipulation that all out-of-plane fluctu- 
ations TOq, = 0; see the Hamiltonian in the next section. 
Therefore, before proceeding to get critical anisotropics. 
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the in-plane vortex structure, as defined using the ring 
variables, was relaxed to a local minimum energy config- 
uration, consistent with equation This was achieved 
by an iteration procedure in which each ring angle 
was adjusted to point along the effective field of its neigh- 
boring rings, scanning through all the active rings, until 
the magnetization and energy changes fell below a de- 
sired precision. Typically, such a relaxation procedure 
only makes very minor changes (of the order of 1% or 
less) in the spin components at sites close to the vortex 
core. Farther from the core, the relaxed spin angles are 
well-described by the continuum formula, Eq. ||2Jl. 



3. Zero-Frequency Dynamic Mode 

Based on the definitions of the rings, it is straightfor- 
ward to determine which pairs of rings are neighbors, 
denoting them as {a, a'). Assuming a symmetrical dy- 
namical mode of oscillation, where the spin fluctuations 
are the same for all sites on a ring (FM model), the total 
system Hamiltonian is written H = Hint + Hgeii, where 
the ring interaction and seZ/ energies are 



^^int — JS ^ ^ Cqq' 



Xr 



{a.a') 



l-ml,cos{<i>l-^l,) , (7) 



(8) 



The sum in Q is only over the neighboring pairs of rings, 
and the sum in ||SJ) is over all rings. 

The classical dynamics for the in-plane spin fluctua- 
tions follows from^ 



dH 
drua 



(9) 



with a similar equation for the conjugate coordinate m^. 
This equation is sufficient to determine when a zero- 
frequency mode occurs, locating a critical anisotropy. 
The equation can be written in a matrix form. 



(10) 



The mode frequency goes to zero when the determinant 
of J- goes to zero, which gives the condition for finding a 
critical anisotropy. This was performed by setting up the 
matrix J- numerically for a number of rings N, and using 
a secant search method, adjusting A or d while holding 
the other fixed, until a zero determinant resulted. 



4- Critical Anisotropies-Numerical Results 

For a total number of rings < 3, one can easily 
write down explicitly the total ring energy and matrix 



TABLE I: Hexagonal lattice critical anisotropy parameters 
for a relaxed vortex on a vacancy as a function of the number 
of rings N and system radius R. 

N R/a Xcv d 

2 1.732 1.5636211762 -0.38565740698 

3 2.000 1.2077648210 -0.20076947906 

4 2.646 1.1101053643 -0.11918190977 

10 4.359 1.0023063573 -0.00297428015 

11 4.583 0.9962392581 0.00492191196 
20 6.245 0.9800726402 0.02721911999 
100 15.000 0.9738836510 0.03669965239 
400 30.643 0.9738607689 0.03674491810 
1500 60.108 0.9738607680 0.03674492096 
3000 85.497 0.9738607680 0.03674492096 



J^, and determine critical anisotropics analytically, as was 
done by Zaspel et && for the square lattice vacancy prob- 
lem. The vacancy-influenced critical exchange anisotropy 
6cv = 1 — Act, is found for null single- ion anisotropy, d = 0. 
The vacancy-influenced critical single-ion anisotropy dcv 
is found for null exchange anisotropy, 5 = or equiv- 
alently, A = 1. The results from a C-program for the 
A^-ring calculation agreed with the analytic calculations 
for A^ < 3 on all three lattices, verifying the reliability 
of the program. Obviously, the critical values depend on 
A^, therefore, results were obtained for a range of A^ up 
to the order of 1000 rings until convergence of Ac and 
dcv to 10 significant figures. A smaller number of rings 
(A^ « 400) is necessary to achieve this for the triangular 
lattice, probably because it is denser. A summary of the 
results is shown in Tables |nl and IIIII On the square 
and hexagonal lattices, the calculation gives no solution 
for either critical anisotropy at A^ = 1. On all the three 
lattices, a minimum number of rings is necessary before 
the critical anisotropics even move into the easy-plane 
range (A^ < 1 and dcv > 0). 

The converged critical exchange anisotropics are Scv ~ 
0.0261392320, 0.0455022615, 0.0606587833, for hexago- 
nal, square and triangular lattices, and the corre- 
sponding single-ion critical anisotropics are dcv ~ 
0.03674492096,0.0824221891,0.1608450953. Just as 
found for a vortex centered in a plaquette, the critical 
anisotropics increase with the coordination number, and 
the single-ion critical value is larger than the exchange 
value. The results confirm and expand upon those found 
by Zaspel et al.^ Most importantly, a vortex on a va- 
cancy will tend to be stabilized in the planar form, even 
at fairly weak anisotropy that would otherwise produce 
stable out-of-plane vortices far from any vacancies. The 
numerical values for the square lattice are slightly differ- 
ent here due to using the energetically relaxed in-plane 
vortex profile (not Eq. IS}, and, calculations to a large 
number of rings until convergence, rather than an ex- 
trapolation procedure from a small number of rings. 
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TABLE II; Square lattice critical anisotropy parameters 
for a relaxed vortex on a vacancy as a function of the number 
of rings TV and system radius R. 



N 


R/a 




dcv 


2 


1.414 


1.4137593264 


-0.4065927450 


3 


2.000 


1.2807409281 


-0.3077361144 


4 


2.236 


1.0608891919 


-0.0871050739 


6 


3.000 


1.0227216056 


-0.0346427396 


7 


3.162 


0.9944246293 


0.0089633765 


10 


4.123 


0.9722646499 


0.0469741982 


100 


15.00 


0.9544993839 


0.0824157890 


400 


30.89 


0.9544977384 


0.0824221893 


800 


44.05 


0.9544977385 


0.0824221891 


1500 


60.80 


0.9544977385 


0.0824221891 



TABLE III: Triangular lattice critical anisotropy param- 
eters for a relaxed vortex on a vacancy as a function of the 
number of rings A'^ and system radius R. 



N 


R/a 




dcv 


1 


1.000 


1.8660254038 


-0.8660254038 


2 


1.732 


1.2709004131 


-0.4284919249 


3 


2.000 


1.0669263945 


-0.1342143055 


4 


2.646 


0.9969826115 


0.0068108464 


10 


4.583 


0.9454955504 


0.1397956043 


200 


24.43 


0.9393412168 


0.1608450949 


400 


35.04 


0.9393412167 


0.1608450953 


800 


50.21 


0.9393412167 


0.1608450953 



B. Relaxed Vortices' Energy and Magnetization 

It is interesting to confirm the locations of the crit- 
ical anisotropics by analyzing the total energy and to- 
tal out-of-plane magnetic moment of a relaxed vortex, 
both as functions of anisotropy. A vortex is considered 
as above, at the center of a circular system centered on a 
vacancy. The in-plane and out-of-plane spin components 
both were allowed to vary, until a local minimum energy 
configuration for H was obtained. In actual practice, the 
minimizing configuration was found by using the origi- 
nal lattice spin fields 6*11 and iteratively repointing each 
along the effective local field Fn due to its neighbors. 



Fn = 



dH 



= J 



1+a 



s s. 



2d Si 



(11) 

Scanning linearly through the lattice, each site was up- 
dated in sequence, being reset along the net field due 
partly to some unchanged neighbors and some that have 
already been repointed. This gives slightly faster conver- 
gence than a synchronized global update. 

Due to the planar symmetry of H , if all spins are ini- 
tially set in the xy plane, they will remain precisely in the 
xy plane under this algorithm, even in a situation where. 



for example, a planar vortex would be unstable. There- 
fore, in order to allow for the possibility to relax either 
to a planar or out-of-plane vortex form, nonzero initial 
values — were set, together with initial in-plane 
vortex angles given by Eq. Q. In situations where an 
in-plane vortex is stable, the S*^ spin components decay 
away to zero; conversely, at anisotropy where an out-of- 
plane vortex is stable, the components grow into the 
appropriate out-of-plane vortex profile. The iterations 
were allowed to proceed until the average spin changes 
per site fell below a desired level, typically on the order 
of 1 part in 10^^. By the time the individual spin changes 
reached this size, the energy converged to an even higher 
precision. It is interesting also to note that this iteration 
converges all the more slowly as the anisotropy comes 
closer to a critical value, a characteristic feature related 
to the slow dynamics of the soft mode responsible for the 
in-plane to out-of-plane instability. 

In a finite system, the boundary can either enhance or 
diminish the ability of nearby spins to tilt out of the easy 
plane, affecting the calculated critical anisotropy values, 
as seen above. For all three lattices studied, fairly good 
convergence of critical values took place below a system 
radius R — 50a. Therefore we calculated relaxed vortex 
states for systems with R — 50a; the spins at the edge of 
the system have no constraint from outside, that is, a free 
boundary condition holds. The energy was calculated rel- 
ative to the ground state energy for spins aligned within 
the xy plane, an amount of —JS'^ per exchange bond. In 
addition to the energy, the total magnetic moment was 
calculated, having only a z component: M = X^n'^n- 
Its deviation from zero gives a distinct signature for the 
transition to an out-of-plane state. 

Results are presented here using exchange anisotropy 
only, setting the single-ion anisotropy to zero. The vortex 
magnetic moments M found as a function of A are shown 
in Fig.^ The critical anisotropics obtained by numerical 
relaxation of the vortex structure, and observed in this 
graph, confirm to several digits the critical anisotropics 
determined in Sec. IIII A 41 The same data is replotted in 
the inset of Fig. 01 versus scaled anisotropy A — Acu, using 
the precise critical values from Sec. IIII A4l One sees that 
once A passes Ac«, the magnetic moment grows initially 
proportional to V A — \cv , and fastest on the more open 
hexagonal lattice, and slowest on the denser triangular 
lattice. There is a corresponding extremely soft change 
in the total system energies (relative to the planar vor- 
tex energy Ep) as shown in Fig. [S] verifying that the 
out-of-plane transition is towards a vortex of lower en- 
ergy, once the anisotropy becomes weaker than the crit- 
ical strength. Some typical out-of-plane vortex profiles 
are exhibited in Fig. indicating not only the growth of 
total M with deviation from critical anisotropy, but also 
an associated increase in the vortex radial length scale as 
the anisotropy strength diminishes. Figure |5| also lends 
support for Zaspel et al's explanation for the lowering of 
the critical anisotropy strengths due to the presence of 
the vacancy. The missing bonds reduce the total energy 
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of both types of vortices, however, the in-plane vortex en- 
ergy is reduced more than the typical out-of-plane vor- 
tex energy (except at rather weak anisotropy) because 
the (removed) out-of-plane tilted exchange bonds were 
relatively smaller energy contributions than if they had 
been lying purely in-plane. This is true on all the lat- 
tices studied, the differences between them being caused 
by the different densities (sites per unit area). 



IV. VORTEX- VACANCY INTERACTION 
POTENTIAL 

The above analysis requires that a vortex is actually 
attracted energetically to a vacancy. This is a reasonable 
supposition; the missing bonds clearly reduce the total 
system energy, and especially the exchange energy due 
to the in-plane components if the vortex is centered on 
the vacancy, because the spins near a vortex core are the 
most strongly misaligned ones. A lesser reduction in in- 
plane exchange energy results if the vortex is near the 
vacancy but not exactly centered on it. Therefore one 
can conclude that there should be an attractive potential 
pulling a vortex into a vacancy site. (M61 et al^^ came 
to a different conclusion, estimating a repulsive potential 
between a vortex and vacancy in a continuum model.) 
Due to its intrinsic interest as well as the implications 
for the acceptable interpretation of the vacancy modi- 
fied critical anisotropics, an analysis of the interaction 
potential between a vacancy and a vortex (in a lattice 
system) is called for. Therefore we investigate further 
the vortex-vacancy energetics, especially as a function of 
their separation. 



A. Vortex-on- Vacancy Pinning Energy 

As a first step, which suggests clearly that the potential 
should be attractive, a comparison is made of the ener- 
gies of a single vortex at the center of a circular system of 
radius R, in two obvious cases: 1) the lattice is uniform 
(no vacancy) and the vortex is centered in a plaquette; 
2) the vortex is centered on a vacancy at the system cen- 
ter. By using the energy minimizing iterative method 
[Eq. Hll|l ] just described, the vortex configurations for 
the two cases are easily obtained, and their energies can 
be compared. At the same radii i?, the two defined sys- 
tems do not necessarily have the same number of lattice 
sites or bonds (their lattices are shifted relative to one 
another and the circular boundary), so the comparison 
may require careful interpretation. We proceeded by pro- 
ducing a set of energies as a function of system size and 
compare the two results versus R. In Fig. [7|the compar- 
ison is made for in-plane vortices on systems with A = 
and d = 0; the vortex-on- vacancy energy lies below the 
vortex-in-plaquette energy for all radii. The difference is 
the vortex-on-vacancy binding energy or pinning energy, 
AE. For example, for in-plane vortices on the square 



lattice, AEsqr « 3.178JS^, more or less independently 
of the system size, and independent of the anisotropy 
as long as the vortices are both of the in-plane type. 
The in-plane vortices on hexagonal and triangular lat- 
tices indicate binding energies of Aii^hex ~ 1.937J5^ and 
AEui w 5.174JS'2, respectively (for X ^ d = 0). The 
conclusion is obvious, that if observed in an extended 
system (aside from possibly strong boundary effects), an 
in-plane vortex will lower its energy by ocurring centered 
on a vacancy, rather than in a plaquette. 

In Fig. |S1 a typical comparison is made for out-of-plane 
vortices, using A = 0.99 (greater then Aci, on all lattices 
studied) and d = 0.0; again the vortex-on-vacancy energy 
is lower for all system sizes. (The calculation required a 
larger system than for A = in order to produce sta- 
ble out-of-plane vortex profiles in a finite radius.) For 
the square lattice, the vortex binding energy is reduced 
drastically, to about 0.232 JS"^; this value is expected to 
be a function of the anisotropy parameters. One could 
also consider an alternative situation where the vortex is 
of out-of-plane form when centered in a plaquette, but 
destabilizes to in-plane form if centered on a vacancy 
(whenever Ac < A < Xcv , such as any of the three lattices 
studied with A — 0.93 and d = 0). In this latter exam- 
ple, the vortex-on-vacancy energy is again consistently 
lower than the vortex-in-plaquette energy at the same 
radius, indicating once again an energetic attraction to a 
vacancy. In this case the pinning process also will be as- 
sociated with a significant signature: the complete elim- 
ination of the out-of-plane magnetization component. 

Binding energies for vortex-on-vacancy pinning are 
given in Table IIVI for the three anisotropy parameters 
A = 0, 0.93, and 0.99, corresponding to permanently in- 
plane, out-of-plane transforming to in-plane, and perma- 
nently out-of-plane vortices. The decrease (AM) in to- 
tal out-of-plane magnetic moment as a vortex is pinned 
on a vacancy is also summarized there. {AM is the 
difference of the vortex-in-plaquette and the vortex-on- 
vacancy magnetic moments.) Pinning energies and mag- 
netic moment changes were calculated for systems with 
radius R — 100a, for the full range of easy-plane exchange 
anisotropy, as displayed in Fig. |^ (As explained above, 
at this radius, the results do not have any dependence on 
system size.) At weak anisotropy (A ^ 1 ) the pinning 
energy increases almost linearly with anisotropy strength. 
The pinning energy saturates at its largest value once 
A < Ac for the lattice under study; the pinning energy is 
always highest for the more dense lattices. The magne- 
tization change on pinning is even more interesting. It 
starts at a nonzero value for the isotropic limit (A — 1), 
reaching a sharp maximum exactly at Xcv for the lattice 
under study (transition between free out-of-plane vortex 
and pinned in-plane vortex). Finally AM becomes zero 
at Ac, where the pinning involves only planar vortices. 

The region where the pinning leaves the vortices in 
out-of-plane form (A > Act;) is of special interest; the 
pinning is associated with relatively small energy change 
in combination with a significant magnetization change. 
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TABLE IV: Some vortex-on-vacancy binding energies (in 
units of JS^) and magnetization reductions (in units of S) es- 
timated by comparing to the vortices-in-plaquette structures 
(d = 0), using circular systems of radius R < 300a. 

A A^hcx AMhex Agsqr AMsqr Agtri AMtri 

0.0 1.937 0.0 3.178 0.0 5.174 0.0 
0.93 1.486 21.31 1.807 28.05 2.45 32.42 
0.99 0.224 24.63 0.232 17.27 0.310 14.85 

Such a vortex pinning/depinning transition might be ini- 
tiated due to thermal fluctuations under appropriate con- 
ditions, associated with large out-of-plane magnetization 
fluctuations. One the other hand, close to the planar 
limit (A = 0), the large pinning energy means that vor- 
tices formed thermally will tend to pin on any vacancies 
present. Free vortex density could be greatly reduced if 
there is high enough vacancy density. 

B. Separated Vortex— Vacancy Potentials 

To proceed further, it is necessary to study the energy 
of a vortex separated from a vacancy by some distance 
Tw The energy minimization scheme of Sec. liii 51 is to be 
applied to different vortex-vacancy arrangements. The 
potential is considered here for square, hexagonal and 
triangular lattices. Generally, the interaction potential 
between the two could depend on the direction relative 
to the lattice, as well as the radial separation. 

In order to investigate this, again a circular system 
of radius R is used, with a vortex placed at its center, 
which is the origin of coordinates, (0, 0). The vacancy is 
placed at some position X = {Xx,Xy) which is allowed 
to vary. Alternatively, if the vortex were placed at differ- 
ent positions away from the system center, a boundary 
energy that changes signiflcantly with the vortex position 
would result, regardless of the presence or absence of a 
vacancy. To avoid this complication, it is much simpler 
to fix the vortex position at the system center, and move 
the vacancy around, measuring the total system energy. 
It is not expected that there is a strong energy change 
due to the vacancy approaching the system boundary. 
Therefore, one expects the resulting total energy to give 
a good indication of the intrinsic vortex-vacancy inter- 
action potential. 

What this means in practice is that, with the vortex 
fixed at the origin (center of circle of radius R), all the 
lattice sites are allowed to shift as the vacancy takes on 
a range of positions. This implies that lattice sites peri- 
odically pass by the vortex as the vacancy is moved to 
a sequence of positions. To avoid an undefined in-plane 
spin direction, a lattice site should not be allowed to fall 
squarely on the vortex position; as a result, some choices 
of vacancy position are prohibited. 

For this numerical algorithm, it is simple to specify 
the vacancy position, which influences the positions of 



the rest of the lattice sites. On the other hand, an exact 
specification of a desired vortex position is not possible, 
because the presence of an off-center vacancy near a vor- 
tex skews the vortex spin field, leading to an ambiguity 
in the position of the "vortex center." Indeed, without a 
sufficient constraint, if an initial nearly in-plane vortex is 
placed near a vacancy, and then the relaxation procedure 
already described fSec. lIIIB)l is applied, it is possible that 
the vortex simply shifts onto the vacant site. Therefore, 
a method is needed to at least partially enforce some de- 
sired vortex position. 

The following simple but somewhat arbitrary scheme 
was applied to enforce a desired vortex core position. 
Based on the arrangement of the lattice sites near the 
vortex, which is not necessarily symmetrical, a few spins 
nearest to the vortex are considered as "core" sites. The 
in-plane angles of these core sites are not changed by the 
energy relaxation scheme using Eq. (|ll|l . Their are 
held fixed at the original angles given according to the 
in-plane profile, Eq. however, their out-of-plane spin 
components are allowed to change. The number of core 
sites is taken to be ric = 6, 4 or 3 for hexagonal, square 
or triangular lattices, respectively. While it is somewhat 
artificial and arbitrary to constrain some spins, it was 
found to be necessary to maintain the desired vortex po- 
sition, especially at vortex to vacancy distances less than 
a few lattice constants. This restriction of some spin 
movement may affect the final energy, therefore one can 
only consider this as an estimate of the vortex-vacancy 
potential. However, at a position of high symmetry, such 
as the vortex on top of the vacancy, the constraint has 
no effect. 

Obviously this is a rich problem with a wide choice of 
parameters. One expects the potential to be influenced 
strongly by the choice of anisotropy constants. Here only 
systems with exchange anisotropy are studied; one ex- 
pects the results for single-ion anisotropy to be compara- 
ble. Furthermore, there could be notably different behav- 
iors depending on where the chosen anisotropy constant A 
lies relative to the vortex-on-plaquette critical anisotropy 
Ac and the vortex-on-vacancy critical value Xcv A par- 
ticularly interesting and complex case is when A lies be- 
tween Ac and A^;; the vortex will take an out-of-plane 
form far from the vacancy, and transform into a planar 
vortex as it is set closer to the vacancy. For A smaller 
than both Ac and A^;, the vortex always remains planar. 
Conversely, for A larger than both Ac and Ac«, the vortex 
always remains out-of-plane, although the magnitude of 
its out-of-plane moment will change as it moves onto the 
vacancy position, as already indicated in Table HVl 



1. Vortex Remaining in a Planar State 

For sufficiently strong anisotropy (A < Ac), the vortex 
always remains planar at any vortex-vacancy separation. 
The parameter choice X = d = was used to obtain this. 
A typical result for the energy is shown in the upper- 
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most curve of Fig. ^1 for a square lattice system, with 
radius R = 50a. Curves with essentially the same shapes 
were obtained also on a system with R = 30a, except for 
an overall shift to lower energies due to the usual log- 
arithmic dependence on system size. One expects that 
except for fine details (such as the width or depth of the 
vacancy-vortex potential), the gross features of the po- 
tential interaction should not depend on the choice of the 
lattice (the usual assumption of continuum theory). 

The curves in Fig. 1101 correspond to vacancy positions 
along paths in the (10) direction of the lattice, avoiding 
those points which would cause the vortex to fall ex- 
actly on a lattice site. A fairly deep attractive minimum 
{AE « 3.3 JS*^) is found for the vacancy centered on the 
vortex, with obvious periodic wings due to the lattice. 
Similar results are obtained for vortex-vacancy displace- 
ments along other directions, but with a slightly different 
width of the potential minimum (the vortex-vacancy po- 
tential is not isotropic). Clearly this confirms the energy 
deviation found in Sec. IIV Al and furthermore, gives an 
indication of the radial separation over which the vacancy 
has a strong influence on a vortex (1-2 lattice constants). 
The periodic variations for large vortex-vacancy separa- 
tion demonstrate the relatively weaker local potential ef- 
fects that tend to center a vortex within an individual 
plaquette. If one did not constrain the small set of core 
spins to enforce a desired vortex position off-center in 
a plaquette, then the relaxation would simply produce 
a vortex profile centered in the nearest convenient pla- 
quette, and these periodic variations would be greatly 
reduced. Similar results hold for the hexagonal and tri- 
angular lattices, but with potential depths as found in 
Sec. IIV Al and periodic variations in the wings that are 
distinctive to the lattice. When A < Ac, the total out-of- 
plane magnetization remains at zero for the whole range 
of vortex-vacancy separations. 



2. Vortex Remaining in an Out-of-Plane State 

It was found above that for out-of-plane vortices, the 
pinning energy was much weaker than for in-plane vor- 
tices. Thus it makes sense also to compare the radial 
dependence of the vortex-vacancy potential interaction 
for out-of-plane vortices to that for in-plane vortices. 

As an example of this possibility, consider exchange 
anisotropy constant A = 0.96 for the square lattice, larger 
than the value Xcv ~ 0.9545, which results in stable out- 
of-plane vortices regardless of their position. The po- 
tential curve for separations along the (10) direction of 
the lattice is shown in the lowermost curve of Fig. [TUl 
Only slightly different potential curves result for sepa- 
rations along other directions. The potential curve is 
very smooth compared to that for the in-plane vortices 
(A = 0), which must be due to the relatively low energy 
out-of-plane tilting adjustments of the spins. The corre- 
sponding change in the overall system out-of-plane mag- 
netization is shown in Fig. 1111 As the vortex is moved 



closer to the vacancy, its out-of-plane magnetization di- 
minishes by AM « 305 at the same time that its en- 
ergy is mildly reduced by AE « l.OJS"^, with an overall 
attraction. Of course, the depth of the attraction will 
depend on the particular choice of anisotropy strength. 

If this type of situation holds in a magnetic material, 
(large change in M associated with small energy dif- 
ference) there could be large out-of-plane magnetization 
fiuctuations expected as vortices are alternately pinned 
and freed by vacancies, due to thermal fluctuations. The 
degree of this effect will be greatest for an anisotropy 
parameter A approaching the isotropic limit, A « 1, cor- 
responding to the largest magnetization change together 
with the smallest energy difference between the free and 
pinned vortex states. The two different states, however, 
could be difficult to distinguish in a real medium, because 
both have nonzero (and large) M. 



3. Vortex Transforming From Out-of-Plane to In-Plane 

Alternatively, consider exchange anisotropy A = 0.90 
for the square lattice, which is above Ac ~ 0.7034 and 
below AcD ~ 0.9545. The vortex-vacancy potential is 
shown as the middle curve of Fig. EH At this anisotropy 
strength, the magnetic moment for the separated out-of- 
plane vortex is about M 195'. As the vortex is pulled 
closer to the vacancy, it minimizes its energy by reduc- 
ing the out-of-plane magnetization to zero, transform- 
ing dramatically from out-of-plane form to in-plane form 
fFig. lllfl . In the process, one sees an intermediate overall 
energy reduction of about 2.3J5^. If such a transition 
between the pinned planar and free out-of-plane states 
could be controlled externally by some applied field, such 
a pair of states could have practical applications for data 
storage, especially because of the transition to an easily 
distinguished zero magnetization state. 



V. CONCLUSIONS 

A vacancy has been found to have substantial influ- 
ence on the structure and stability properties of a vortex 
in an easy-plane magnet. The quasi-static analysis pre- 
sented here is valid for FM and AFM models on square 
and hexagonal lattices, and only for the FM model on 
a triangular lattice. Consistent with the stability results 
found by Zaspel et alji when a vortex is energetically re- 
laxed centered on a vacancy, relatively weaker easy-plane 
anisotropy will stabilize it in the planar form, compared 
to that necessary when the vortex is centered in a pla- 
quette far from any vacancy. This result holds as well on 
hexagonal and triangular lattices; the critical anisotropy 
strengths 5cv and dcv are always reduced for vortex-on- 
vacancy compared to those for vortex-in-plaquette, 5c 
and dc- All the critical anisotropy strengths are small- 
est on the hexagonal lattice and largest on the triangular 
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lattice, which suggests that the coordination number and 
lattice density control the effect. 

Using a simple scheme to enforce a desired vortex posi- 
tion, energetically relaxed vortex-vacancy configurations 
were produced. The potential interaction between the 
two has been estimated for the different lattices and 
anisotropy regimes. In all cases studied, the vortex- 
vacancy interaction is attractive. This is important be- 
cause it establishes the vortex-on- vacancy states as legit- 
imate local minimum energy configurations. The depth 
of the potential holding a vortex on a vacancy (pinning 
energy) has been estimated. It tends to be strongest on 
the triangular lattice and weakest on the hexagonal lat- 
tice for in-plane vortices. For out-of-plane vortices, the 
typical depth of the vortex pinning potential tends to be 
much weaker than for in-plane vortices. 

An interesting effect is expected for a system with 
anisotropy strength intermediate between the vortex-in- 
plaquette and vortex-on- vacancy critical anisotropics (for 
example, Ac < A < Xcv)- In the transformation from out- 
of-plane to in-plane type as a vortex moves from free to 
pinned on a vacancy, the out-of-plane magnetization is 
completely eliminated. This could lead to strong magne- 
tization fluctuations, and also be considered as an addi- 
tional signature of vortices and a distinct signature of the 
pinning process. Additionally, the magnetic fluctuations 
associated with pinning and de- pinning of out-of-plane 
vortices may be even stronger: relatively large out-of- 
plane magnetization changes will occur at very low en- 
ergy difference between the pinned and free vortex states. 
The largest magnetization changes in the pinning pro- 
cess occur at the vortex-on-vacancy critical anisotropy, 

A ~ Xcv 

Alternatively, one could consider application of a weak 
out-of-plane magnetic field, which would bias the out- 
of-plane polarizations of generated vortices along that 
preferred direction (generating the "light" cone-state 



vorticesiSiii). An initial low-temperature state of pla- 
nar vortices pinned on vacancies could produce de-pinned 
out-of-plane vortices preferentially with the same polar- 
ization, leading to a large macroscopic magnetic moment. 
This could be a intriguing effect if the de-pinning and pin- 
ning of vortices on vacancies could be controlled at will, 
leading to interesting device possibilities. 

Finally, one can speculate on how an attractive po- 
tential between vacancies and vortices should affect the 
statics and dynamics of the BKT topological transition. 
The MC results by Leonel et ali^ indicated that a small 
percentage of spin vacancies lowers the transition temper- 
ature, which they interpreted to result from a repulsive 
vacancy-vortex potential (under an assumption of a cer- 
tain vortex deformation caused by a vacancy). The ex- 
planation requires that the repulsion produces a greater 
free- vortex density at a lower temperature, thereby lead- 
ing to more disorder and a lower transition temperature. 
Actually, we found here that the vacancies are likely to 
attract and pin any vortices/ antivortices that are form- 
ing, with fairly strong pinning energies for the planar 
model, and associated lower vortex formation energies. 
An improved continuum model for vacancy-vortex in- 
teraction, without the assumption of a global vacancy- 
induced vortex deformation, has verified this attractive 
potential. "'^^ At temperatures well below the energy scale 
of JS^, nearly all formed vortices should be pinned. Even 
so, the pinned vortices will lead to long range disorder 
and the BKT transition, but at a lower temperature, es- 
pecially because their formation energy is much lower 
than for unpinned vortices. Further detailed MC studies 
could be interesting if they could demonstrate the respec- 
tive roles of the free versus pinned vortices, by analyz- 
ing their number densities. This question may also be 
answered by studying the dynamic correlations: the dy- 
namic signatures of free vortices should be distinct from 
those due to pinned vortices. 
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FIG. 1: Vacancy (x) on a square lattice, and surrounding 
sites which define the (i, j) coordinates of the rings (solid 
points in first octant arc the principal sites). Numbers indi- 
cate the symmetry related sites of some of the first few rings. 
Rings with principal site on one of the symmetry lines at 0° 
or 45° have four sites; all others have eight. 
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FIG. 2: Vacancy (x) on a hexagonal lattice and definitions 
of the (j, j) coordinates of some of the first few rings. Rings 
with principal site on one of the symmetry lines at 0° or 60° 
have three sites; all others have six. 
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FIG. 3: Vacancy (x) on a triangular lattice and definitions 
of the (j, j) coordinates of some of the first few rings. Rings 
with principal site on one of the symmetry linos at 0° or 30° 
have six sites; all others have twelve. The dashed lines exhibit 
some of the intra-ring or self-exchange bonds. 




FIG. 4: Total out-of-plane magnetization of a numerically re- 
laxed vortcx-on- vacancy versus exchange anisotropy constant 
A, for hexagonal, square and triangular lattice systems of ra- 
dius -R = 50a with free boundaries. The inset shows the same 
data versus shifted exchange anisotropy constants A— A 
the critical points. 
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FIG. 5: Energy of a numerically relaxed vortcx-on-vacancy 
relative to in-plane vortex-on- vacancy energies i5vv(0) versus 
shifted exchange anisotropy constant A — \cv, for hexagonal 
{E^^{Q) = 7.143JS2), square {E^^{0) = 13.15JS^) and trian- 
gular (Svv(O) = 23.46 JyS^) lattice systems of radius R = 50o 
with free boundaries. 




FIG. 6: Radial dependences of the out-of-plane component 
of numerically relaxed vortex-on-vacancy structures, at A = 
0.975, 0.980, 0.985, 0.990, 0.995, on a hexagonal lattice system 
of radius R = 50a with free boundaries. 
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FIG. 7: Comparison of vortex-in-plaquette energies (dotted) 
with vortex-on-vacancy energies (solid) for in-plane vortices 
(A = d = 0) on circular systems of radius /?, showing a sub- 
stantial pinning energy. 
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FIG. 8: Comparison of vortex-in-plaquette energies (dot- 
ted) with vortex-on-vacancy energies (solid) for out-of-plane 
vortices (A = 0.99, d = 0) on circular systems of radius _R, 
showing a very weak pinning energy (See Table HVI '). 
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FIG. 9: The changes in vortex out-of-plane magnetie moment 
(AM) and in vortex energy {AE, in inset) when a vortex is 
pinned on a vacancy (changing from centered in a plaquette to 
centered on a vacancy) , as functions of easy-plane anisotropy 
parameter on hexagonal, square and triangular lattices. 




FIG. 10: Vortex-vacancy total energy as a function of their 
separation rvv, calculated on a square lattice system of ra- 
dius R — 50a for indicated exchange anisotropics (A = 0, 
in-plane; A = 0.90, transition from in-plane to out-of-plane 
with A = 0.96, out-of-plane vortices). With 

the vortex at the origin (0,0), the vacancy was placed at a 
sequence of positions in the (10) direction. 
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FIG. 11: Vortex-vacancy total out-of-plane magnetization 
as a function of separation Tw, for square lattice systems of 
radius R = 50a as in Fig. [HIl (M = for all rw at A = 0.) 



